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Abstract 

We present the results of an extension of our previous work on large-scale simulations of 
dynamically triangulated toroidal random surfaces embedded in with extrinsic curva- 
ture. We find that the extrinsic-curvature specific heat peak ceases to grow on lattices with 
more than 576 nodes and that the location of the peak Ac also stabilizes. The evidence 
for a true crumpling transition is still weak. If we assume it exists we can say that the 
finite-size scaling exponent ^ is very close to zero or negative. On the other hand our 
new data does rule out the observed peak as being a finite-size artifact of the persistence 
length becoming comparable to the extent of the lattice. 
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The theory of 2d fluid random surfaces embedded in R^, with an ex- 
trinsic curvature term (bending rigidity) in the action, has received consid- 
erable analytical and numerical attention in the last decade [1, 2, 3]. In 
[3] we presented the results of a large-scale Monte Carlo simulation of a 
dynamically-triangulated torus in with up to 576 nodes, corresponding 
to 1152 triangles. Although we observed a rapid crossover from a crumpled 
regime for A < Ac to a smooth regime for A > Ac, where A is the extrinsic 
curvature coupling constant and Ac ~ 1.425, it was not at all clear whether a 
true continuous thermodynamic phase transition separated the two regimes. 
In fact several alternative interpretations of the data were discussed in [3]. 
Perhaps the simplest possibility, advocated in [4], is that the persistence 
length ^ describing the exponential decay of the normal-normal two-point 
function in the crumpled (disordered) regime simply reaches the finite size 
of the system at Ac. In this case the observed smooth regime would be a 
finite-size artifact with the true continuum theory really being crumpled for 
all couplings A, in accordance with perturbative analytical results [5, 6, 7, 8]. 
Since ^ grows exponentially with A, according to the one-loop beta-function, 
this interpretation would imply that Ac diverges logarithmically with system 
size A'^. To resolve this issue and to gain further insight into the model it 
was clearly desirable to extend the numerical simulations to larger lattice 
sizes and to clarify the influence of flnitc-size effects. In this short letter we 
present an extension of our previous work to include toroidal lattices 1152 
and 2304 nodes. 

As in [3] we study the theory defined by the action 

S = Saauss + XSe=J2 C^J{Xt - X^y + A^(l - • nf) , (1) 

where Cij is the adjacency matrix, Xf is the position in R^ of node i 
{i — 1,..,N) and is the normal vector to a triangle Ic in the cellular 
decomposition of a lattice discretization of a torus. The discretization 5'^ of 
a continuum extrinsic curvature term takes support on the edges (links) of 
the lattice and is known as the (discrete) edge extrinsic-curvature. The sim- 
ulation consists of a standard Metropolis algorithm for the updating of the 
nodes and a DTRS- algorithm [9, 10] to sweep through the space of trian- 
gulations. The basic flip move is attempted on randomly chosen links. After 
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a set of 3A^ flips are performed, 3A^ randomly selected embedding coordinates 
are updated by random shifts from a flat distribution. 

The observable of most direct physical interest is the edge extrinsic- 
curvature specific heat 

C{X) ^^{<Sl>-< Se >') . (2) 

This exhibits a peak at a coupling Ac which depends on the exact discrete 
form of the action chosen [11, 12, 13, 14, 15, 16, 2, 3]. In [3] we found 
that the maximum value of the specific heat grows with the system size as 
Cmax = AN'^, with uj = 0.06 ± 0.05. In our new series of simulations on 
lattices with 1152 nodes we ran 54 million sweeps at A = 1.425, 21 million 
sweeps at A = 1.430 and 18 million sweeps at A = 1.435 ^ On the data 
from these three points we use multi-histogram reconstruction [17, 18]. This 
works well in that three different reconstructions give coherent results. On 
lattices of 2304 nodes we have poorer statistics. We ran 17 million sweeps 
at A = 1.425 plus approximately 5 million sweeps at A = 1.40, 1.42 and 1.43 
as a consistency check. On the 2304 lattice histogramming does not work 
well. This is to be expected since the statistics are not good enough for such 
a large lattice. Still we have checked that our measurements at A = 1.425 
give consistent results, that the error estimate is reliable and that we are, 
with good accuracy, at the peak of the specific heat. In all these simulations 
required the equivalent of approximately one year of CPU time on an HP 
9000 (720 series) workstation. 

The specific heat peak for N — 576, 1152 and 2304 is shown in Fig. 1. 
In Table 1 we give our results for the maximum of the specific heat and the 
associated coupling Ac function of . 

Clearly the maximum of the specific heat curve Cmax is effectively con- 
stant for surfaces with 576 or more nodes. The (pseudo)-critical coupling Ac 
is also constant for N — 576 and above. With the present data we can defi- 
nitely exclude the presence of a divergence in the specific heat. The growth of 
the specific heat peak observed on small lattices [11, 12, 13] does not refiect 

^ Since the autocorrelation time t is of order 400, 000 sweeps on the 1152 lattice these 
runs have at least 45r measurements. 

^We have reanalyzed the data presented in reference ([3]), using a different method of 
weighting relative errors when combining histograms. Thus, some of the errors quoted 
here are smaller than the respective uncertainties in reference ([3]). 
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Ac 


36 


3.484 ± .008 


1.425 ± .035 


72 


4.571 ± .015 


1.410 ± .015 


144 


5.37 ± .08 


1.395 ± .017 


288 


5.55 ± .05 


1.410 ± .015 


576 


5.81 ± .06 


1.425 ± .010 


1152 


5.69 ± .04 


1.425 ± .010 


230 1 


5.75 ± .10 


1.125 ± .010 



Table 1: The maximum of the specific heat and its position, with errors, for 
different lattice sizes. 

true asymptotic behaviour. These results also invalidate the interpretation 
raised in the introduction [4, 3]. A one-loop renormalization group calcu- 
lation shows that the persistence length ^ grows with bending rigidity A as 
exp( J^A). Equating ^ with the spatial extent of the lattice A^^/'^*", where din 
is the intrinsic Hausdorff dimension of the lattice, one sees that they become 
comparable at a coupling Ac ~ In(A^). In the continuum limit A^ oo, 
Ac diverges. Since, for reasonable values of din, we do not see the increase 
in Ac with N predicted by the above relationship we can state with some 
confidence that the origin of the observed specific heat peak is not explained 
by the persistence length becoming comparable to the extent of the lattice. 

In [3] we also discussed other, more subtle, possibilities that could account 
for the observed behaviour of C(A) without invoking a phase transition. One 
was based on the analogy between the present model and the 0(3) sigma- 
model in 2d [19]. This model is also asymptotically free and consequently 
disordered at all non-zero temperatures. Yet numerical simulations show a 
distinct peak in the specific heat which grows for small lattices and then 
saturates, just as we find in the model of a rigid string treated here. This 
may be seen in Fig. 2 where we have plotted the specific heat C{P) for the 
two dimensional 0(3) model. The simulations were done on square lattices 
of volume A^ = 16, 25, 64, 100, 900, 2, 500, 4, 900 and 10, 000 using the Wolff 
algorithm [20]. For each point of the A^ = 25,100 and 10,000 lattices we 
used 100,000 measurements. We took a measurement every time the Wolff 
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clusters updated a volume exceeding 30 times the volume of the lattice ^. 
For each point of the = 16, 64, 900, 2, 500 and 4, 900 lattices we used 
20, 000 measurements. We took a measurement every time the Wolff clusters 
updated a volume exceeding 3 times the volume of the lattice. It is very clear 
that the peak levels off quickly for > 100 and that "/?c" is not increasing 
with the size of the lattice. Measurements of the asymptotic value of C(/9) 
have been reported in the past [21, 22]. The authors of [23, 24, 25] explain 
the peak as the excitation of an extra degree of freedom, the so-called a- 
particle [26]. The would-be transition occurs when the mass of the a-particle 
becomes comparable to the inverse correlation length of the 0(3) model. It 
may be that there is a similar interpretation of the observed peak of C(A). 

We are currently histogramming our data to examine as well the be- 
haviour of the complex zeroes of the partition function when A is allowed 
to become complex. For SU(2) lattice gauge theory, which also exhibits a 
specific heat peak without an associated phase transition, it has been shown 
that there are complex zeroes which are near the real axis but do not con- 
verge to it in the infinite- volume limit [17]. High-temperature expansions 
also indicate that the 0(3) model susceptibility has a complex singularity 
near the real axis [27] . 

It is still possible that there is a true continuous phase transition, the 
crumpling transition, occurring at Ac. Assuming a continuous transition, a 
standard finite-size-scaling argument only tells us that u = ^ < 0, where 
a is the exponent governing the divergence of the specific heat, 0(A) ~ 
|A — Acp*^, ly is the analogous exponent for the correlation length and d is 
the intrinsic Hausdorff dimension of the surface. In other words there may 
be a cusp singularity at Ac as, for example, in the case of the superfluid (A) 
transition in He^ [28, 29], for which a = -0.0127 ± 0.0026. Since we do 
not have a measurement of vd, which may even be rather large, we have 
no reliable idea of the exponent a itself. Generally speaking one finds that 
second order transitions on fixed lattices become higher order on dynamical 
lattices, as for example in the case of the 2(i-Ising model [30, 31]. Since there 
seems to be a 2nd order crumpling transition for non-self-avoiding tethered 
(fixed-triangulation) surfaces [33, 34, 35, 36], it would be consistent for the 
transition to be higher than 2nd order when the model is coupled to gravity. 

^For the N = 16, 64, 900 and 4, 900 lattices the integrated autocorrelation times were 
between 1 and 2 Wolff updatings of the entire lattice 
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N 


F[K] 


F[q - 6] 


576 


5.71 ± .08 


8.39 ± .04 


1152 


5.59 ± .05 


8.32 ± .03 


2304 


5.70 ± .10 


8.37 ± .06 



Table 2: The mean square fluctuations of the extrinsic Gaussian curvature 
/C and the defect coordination number q — 6, with errors, for different lattice 
sizes. 

All told our work gives only weak evidence for a continuum crumpling 
transition. The strongest evidence in favour of such a transition, at present, 

is the scaling behaviour of the string tension and mass gap reported in [2]. 
This highlights the need for more extensive measurements on these important 
observables. 

We have also measured the fluctuations of the extrinsic Gaussian curva- 
ture |/C|, defined as the average magnitude of the deficit angle at each vertex 
as measured in the embedded space. Likewise, we have also computed the 
fluctuations of the mean defect coordination number |g — 6|, which is propor- 
tional to the intrinsic Gaussian curvature. In reference [3] we had observed 
that fluctuations of these observables were quite large near the coupling Ac 
but then drop quite dramatically for slightly higher A. We find that on larger 
lattices the fluctuations of these observables at Ac also do not grow with N; 
thus their behaviour does not provide unequivocal evidence of the presence 
of a phase transition. In Table 2 we give the mean-square fluctuations of 
both observables. 

Finally we note the behaviour of the gyration radius at Ac (which we 
will take as being 1.425 in the following). For large A (> 2), the scaling of 

2 

R{N) N'^H with gives a Hausdorff dimension close to 2 (as we expect 
for flat surfaces). In the crumpled region the Hausdorff dimension rapidly 
increases with diminishing A. We had pointed out in [3] that finite size 
effects were relevant in the sector close to Ac and that we could not estimate 
a reliable number from the lattice sizes analyzed. Here the largest lattice 
we simulated (A^ = 2304) does not give useful data, since the error in R 
is too large, but on the 1152 and 576 node lattices we get a fairly precise 
estimate of R, which allows us to estimate for the Hausdorff dimension at the 
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pseudo-critical point Ac the value dn = 4.35 ± .3. This is an intriguing result, 
since 4 is the extrinsic Hausdorff dimension of a class of branched polymers, as 
constructed, for instance, in references [37] . Such configurations are expected 
to dominate the string functional integral for large embedding dimension D. 
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Figure Captions 



1 The edge extrinsic-curvature specific heat C(A) as a function of A. 

Multi-histogram reconstructions with errors are shown for = 576 
(long and short-dashed lines) and = 1, 152 (solid lines). Four indi- 
vidual data points are also shown for = 2,304 (solid circles). One 
sees that the specific heat peak has saturated - it is not growing with 
the system size N above 576. 

2 The specific heat C{(3) of the two-dimensional 0(3) non-linear sigma 
model as a function of /? for lattice volumes N — 16, 25, 64, 100, 900, 
2, 500, 4, 900 and 10, 000. The peak saturates quickly for N > 100 and 
"/3c" does not increase with the volume. 
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